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Abstract 

CO ■ We study the basic quantum mechanics for a fully general set of 

^^ Dirac matrices in a curved spacetime by extending Pauli's method. 

We further extend this study to three versions of the Dirac equation: 
IO ■ the standard (Dirac-Fock-Weyl or DFW) equation, and two alterna- 

^— . . tive versions, both of which are based on the recently proposed linear 

tensor representations of the Dirac field (TRD). We begin with the 
qq ■ current conservation: we show that the latter applies to any solution 

^^ . of the Dirac equation, iff the field of Dirac matrices 7^ satisfies a spe- 

cific PDE. This equation is always satisfied for DFW with its restricted 
choice for the 7^ matrices. It similarly restricts the choice of the 7^ 
matrices for TRD. However, this restriction can be achieved. The 
frame dependence of a general Hamiltonian operator is studied. We 
show that in any given reference frame with minor restrictions on the 
spacetime metric, the axioms of quantum mechanics impose a unique 
form for the Hilbert space scalar product. Finally, the condition for 
the general Dirac Hamiltonian operator to be Hermitian is derived in 
a general curved spacetime. For DFW, the validity of this hermiticity 
condition depends on the choice of the 7^ matrices. 

Key 'words: Dirac equation, gravitation, current conservation, Her- 
mitian Hamiltonian, tensor representation 



1 Introduction 

1.1 Aim of this work 

Quantum-mechanical effects of the classical gravitational field are currently 
the clearest consequence of the interplay between gravitation and the quan- 
tum. The effects which have been observed [151 1221 123 E7J HI] had been 
previously predicted by using the non-relativistic Schrodinger equation in 
the Newtonian gravity potential [23 ED E], and are still described by this 
same non-relativistic approximation [281 [371 02] . However, one expects that 
the precision will increase so that, in the future, the corrections brought by 
the Dirac equation in a curved spacetime (see e.g. Refs. [TT | [20 | [26 | [29 | li m i2]) 
should become detectable. The standard extension of the Dirac equation to 
a curved spacetime is due to Fock and to Weyl, and will be referred to as 
the "Dirac-Fock-Weyl" (DFW) equation. In addition to this, two alternative 
versions of the Dirac equation in a curved spacetime were tentatively pro- 
posed in a recent work [3J. While waiting for an increase in the experimental 
accuracy, one would like to determine the predictions of the Dirac equation 
for quantum mechanics in a curved spacetime, and to check if the different 
versions of it might be experimentally distinguishable. 

A basic feature of the Dirac equation is that its coefficients, the Dirac 
matrices 7 M , have to satisfy the anticommmutation relation corresponding to 
the Lorentzian metric g^ u on the spacetime V: 

yY + 7V = 2<ri4, n,ve{0,...,3} (l 4 = diag(l, 1,1,1)). (1) 

[Here {g^ u ) is the inverse matrix of (g^u)-] Therefore, in a curved spacetime, 
or already in a flat spacetime in general coordinates, the Dirac matrices 
7 M (X) depend on the spacetime point X, as does the metric. It follows that 
there is a continuum of possible choices for the field ■j^(X), all satisfying 
the fundamental anticommutation relation (TjQ). The point- dependence of the 
metric also implies, as we will show, that there is an infinity of a priori equally 
valid candidates for the Hilbert space scalar product. 

Together with the standard equation or Dirac-Fock-Weyl (DFW) equa- 
tion, the two alternative equations [3] provide us with three versions of the 
Dirac equation in a curved spacetime, which are a priori inequivalent. The 
aim of the present work is to study the basic quantum mechanics: current 



conservation, Hilbert space inner product, and hermiticity of the Hamilto- 
nian, for these three gravitational Dirac equations. To reach this goal, in 
particular to study the influence of the possible choices for the field of Dirac 
matrices 7^(X), it is necessary to be able to use any possible choice of the 
latter. This is achieved by using the "hermitizing matrix" A of Bargmann 
[S] and Pauli [531 US] , which allows one to define the current and the scalar 
product for a generic (ordered) set (7^) of Dirac matrices [3j. To our knowl- 
edge, until now, the possible choices for the fields 7 M (X) and the Hilbert 
space scalar product have not been systematically investigated, even for the 
standard (DFW) equation. [For the latter, the different fields 7 M (X) arise 
due to different choices of the orthonormal tetrad field in Eq. (|T5|) below.] 
However, it is not a priori obvious that these choices have no effect. We 
will show that in fact they have an essential effect. First, we will show that 
the axioms of quantum mechanics impose a unique form for the scalar prod- 
uct. The hermiticity condition has not been investigated in a general setting, 
again even for the standard (DFW) equation. The recent work of Leclerc [21] 
studies the hermiticity of the Hamiltonian without explicitly stating restric- 
tions on the coordinate system, but it restricts consideration to special kinds 
of tetrad fields, and the validity of its result is not general: we will show that, 
for DFW, the validity of the hermiticity condition depends on the choice of 
the fields 7 M (X). This indicates the presence of a uniqueness problem. 

For example, this uniqueness problem arises for the DFW equation in a 
Rindler spacetime, which is a Minkowski spacetime using Rindler coordinates 
[30] . In a Rindler spacetime, it is usual to choose a tetrad that is independent 
of the time coordinate. However, if time dependent tetrads are allowed, the 
time evolution of the quantum mechanical states can fail to be unitary. How 
does one exclude the use of the Minkowski tetrad which is time dependent in 
a Rindler spacetime? Furthermore, there is no known method for choosing 
this tetrad that can be generalized to general spacetimes, such that the time 
evolution will be unitary. The approach taken in this paper is complemen- 
tary to using tetrads, and questions such as the existence and uniqueness 
of a Hermitian Hamiltonian operator for a given Dirac equation — which is 
essential for energy eigenfunction expansions in both time-independent and 
time-dependent quantum mechanics — are resolved more directly by consider- 
ing admissible local similarity transformations. (See Section HJ) Clearly, the 
present investigation of the first quantized Dirac theory in a curved space- 
time is a preliminary step to second quantization in a curved spacetime. 



1.2 Status of the two alternative Dirac equations 

The two alternative versions of the Dirac equation in a curved spacetime were 
got by using directly the classical-quantum correspondence [3]. One version 
obeys [U |3] the equivalence principle in the sense which is standard [13] 
in this context. Namely, it automatically coincides with the flat-spacetime 
Dirac equation "in a local freely- falling frame," i.e., in a coordinate system in 
which, at the point X considered, the metric tensor reduces to the standard 
form i]^ u [with matrix (r/^) = diag(l, — 1, — 1, — 1)], and the metric connec- 
tion vanishes. In contrast, for the DFW equation, the equivalence principle 
can only be established for an anholonomic frame and a spin connection [3]. 
The other alternative gravitational Dirac equation has a preferred reference 
frame, although it can be rewritten in a generally-covariant form [3]. Both 
of these alternative versions require the tensor representation of the Dirac 
(TRD) field, i.e., the Dirac theory with a four- vector wave function and with 
the set of the components of the 7 M matrices building a third-order tensor 
[I]. In the TRD theory, the Dirac equation has the same linear form and a 
similar Lagrangian as the DFW equation, \j potentially giving rise to both 
first and second quantized theories. 

In a Minkowski spacetime in both Cartesian and affine coordinates, the 
TRD theory with constant Dirac matrices has been proved [1] to be quantum- 
mechanically fully equivalent to the genuine Dirac theory. See also Subsect. 
1.1 in Ref. [3] for a summary of the argument. Like the solutions of the 
genuine Dirac equation, the single particle TRD solutions have only two spin 
polarizations (up and down) which makes them spin-half wave functions. 
Moreover, just as with the genuine Dirac equation, it is straightforward to 
extend the single particle theory to a canonical second quantized fermion 
theory (see Subsect. 4.3 of Ref. [38]), with the fermion field operator built 
on the normalized single particle and antiparticle TRD solutions. 

A global field of Dirac matrices 7 M satisfying Eq. (JT]) exists (for DFW 
and for TRD as well) if and only if the spacetime admits spinor structure 



1 See Eq. ^ below for the common form of the Dirac equation and see Eq. (21) in 
Ref. [5] for the Lagrangian. 



It thus suggests itself to extend the investigation to curved spacetimes 
that admit spinor structure. We will first present the three Dirac equations 
in a common framework and show that appropriate hermitizing matrices al- 
ways exist (Sect. |2]). Then we shall discuss the definition of the current, the 
condition for its conservation, and the ways to fulfil this condition (Sect. [3]). 
Finally, we shall write the Hamiltonian and emphasize its frame-dependence, 
define the scalar product, and characterize the hermiticity condition for the 
Hamiltonian (Sect. H|). In conclusion we summarize the theorems (proved in 
this paper) which are common to all three Dirac equations (Sect. |5]). 



2 Three Dirac equations in a curved space- 
time 

2.1 The Dirac equation with three different connec- 
tions 

The three gravitational Dirac equations discussed in the llntroductionl have a 
common form: 

7"£> M V = -imif), (2) 

where 7 M = 7 M (X) (n = 0, ..., 3) is a field of 4 x 4 complex matrices defined 
on the spacetime manifold V, satisfying the anticommmutation relation (JTJ); 
and where i = a/— 1, m is the rest-mass of the particle (setting H = 1 = c: 
otherwise, replace m by mc/h), ip is either a quadruplet of four scalar fields 
(DFW) or a four- vector field (TRD), and D^ is a covariant derivative, asso- 
ciated with a specific connection. 

For the two alternative equations based on the tensor representation of the 
Dirac field (TRD), this is an affine connection (i.e., a connection associated 
with the tangent bundle [HUB]): 

(D^y = d,r + a^. (3) 

In the TRD equations, the affine connection is extended to the complexified 
tangent bundle. More precisely, for one of the two TRD equations, henceforth 
denoted TRD-1, this is simply the Levi-Civita connection. That is, the 



A£ 's are the (second- kind) Christoffel symbols associated with the spacetime 
metric g^ v : 

** = &}■ W 

This is the one which obeys the equivalence principle, in the sense stated in 
the [Introduction! For the other TRD equation (TRD-2), the connection A 
is defined from the spatial Levi-Civita connection in an assumed preferred 
reference frame (a congruence of observers [13] , here endowed with a preferred 
time coordinate [U [6] , denoted E. In any coordinates adapted to E, we have 

{0 ifz/ = 0or/i = 0orp = 
(5) 
G\ k if v — j and ft = k and p = I G {1, 2, 3}, 

the Gj k 's (j, k, I G {1, 2, 3}) being the Christoffel symbols of the spatial met- 
ric h [131 [23J |27] in that preferred frame E. In a general coordinate system, 
the Ap 's are obtained from the coefficients (JSJ), by using the transforma- 
tion law of affine connections |17j : it is proved in Ref. [3] that this does 
define a unique connection (which is torsionless) . Equation TRD-2, being 
based on that "preferred-frame connection," has clearly a more speculative 
character from the physical point of view. Note that, in the case of a static 
(and non-flat) spacetime, we do have one preferred reference frame with a 
preferred time coordinate [3]. Also note that, in the case of a flat space- 
time, special relativity does apply. Thus, in that case, we must find that we 
may take any inertial frame as "the preferred frame E." And indeed, if we 
take for E any inertial frame, the connection defined by (jSJ) coincides with 
the Levi-Civita connection associated with the flat metric (with, in particu- 
lar, all A^ 's being zero if the coordinates are Cartesian). In other words, 
the two equations TRD-1 and TRD-2 coincide in the case of a flat spacetime. 

For the DFW equation, the connection is the "spin connection" acting 
on the trivial bundle V x C 4 . g It is built from the "spin matrices" T M [12J: 



2 It is well known that a given spacetime V need not admit a spinor structure. It was 
proved by Geroch that a four-dimensional noncompact spacetime admits a spinor structure 
if and only if it admits a global tetrad field [19]. In that case, both spinor bundles and 
tangent bundles are trivial [TH1 [H] ■ Penrose and Rindler argue that these are the only 
spacetimes of interest [35] . Then it is easy to show that one may define a global held of 
Dirac matrices 7^ satisfying Eq. (JT]) above, for DFW and for TRD as well. 



we shall take the positive sign, 

d^ = d^ + r>, (6) 

so that ([3]) takes the form (^ if we associate matrices T M to connection 
coefficients (or conversely) by 

W'p - A^. (7) 

However, except in the Majorana representation where the T^'s are real, the 
spin connection matrices T^ of the DFW equation are generally complex, and 
depend on the set of fields (7^) (// = 0, ..., 3). 

The transformation of the Dirac equation (|2j) for a coordinate change 
depends also on the version: for the two TRD equations [3J, the wave function 
i[> is a spacetime vector (or four- vector) , hence it transforms thus: 

^ = 14 w = u* u r), L ^^> ( 8 ) 

and the threefold array of the components of the Dirac matrices, 7 PP = (7 M ) P ^, 
builds a (f ) tensor, thus 

^ = ^DfL-\ or 7 r = L^L" T (L- 1 )^ 7 7. (9) 



In contrast, for the DFW equation |12j . the wave function transforms like a 
scalar: 

^{{x ,v )) = r((x u )l (10) 

and the set of the 7^ matrices transforms like a four-vector: 

l ni = L^Y- (11) 

Either of these two transformation modes leaves the Dirac equation fl2]) form 
invariant after a general coordinate change: for the TRD equations, each 
side of (J2j) is a four-vector; whereas, for DFW, it is an object with four scalar 
components. 



2.2 A common tool: the hermitizing matrices 

To define the current, we shall use the hermitizing matrix for a general set 
(7^) of Dirac matrices. This is a nonzero 4x4 complex matrix A such that 

A*=A, (AYY = AY /i = 0,...,3, (12) 

where M< = M* T denotes the Hermitian conjugate of a matrix M. Moreover, 
we shall see that the Hamiltonian operator associated with (J2J) depends on 
the field of matrices a M , with 

a ee //.A a* ee 7 V7^ °. (13) 

Therefore, to define a relevant scalar product, we will also need a hermitizing 
matrix, denoted B, for the set (a M ): 

B f = .B, (St*")* = Ba" /i = 0,...,3, (14) 

and we will need that it be positive- definite. The existence of A and B is 
ensured by the following result [I]: g 

Theorem [4j-6. Fix any point X in the spacetime. For any set of ma- 
trices 7 M satisfying the general anticommutation formula (QJ), there exists a 
hermitizing matrix A for the matrices 7 M . The matrix A is nonsingular and 
unique, up to a real scale factor. Similarly, a nonsingular hermitizing matrix 
B = Aj° for the a^ 7 s exists and is unique, up to a real scale factor. If, 
furthermore, the coordinate system is an admissible one, i.e., if goo > and 
the 3x3 matrix (gjk) (j, k — 1,2, 3) is negative definite, then B = A^ is 
either a positive or negative definite matrix. The sign of the matrix A can be 
chosen such that B = A'y is a positive definite matrix. 

The proof of this theorem in Ref. [1] is directly valid for TRD, because 
that proof uses the tensor transformation of the Dirac matrices, Eq. §§§ here, 
but this is not in an essential way. For DFW, the field of the matrices 7 M 
satisfying the anticommutation relation (pQ) is defined from an orthonormal 
tetrad (u a ), with u a = a\ £- (a = 0, ..., 3), by P2 122 128] 

I 11 = a\ ^ a , (15) 



3 The existence of a hermitizing matrix A had been already proved by Pauli [33l [34] , 
though in a less general case and with less complete results (see Ref. [4]). 
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where (7^*) could be any set of constant Dirac matrices, i.e., obeying Eq. (pQ) 
with (g^) = {ji^y)^ 1 = (t)^) = diag(l, — 1, — 1, —1). It is immediate to check 
that, due to the fact that the matrix a = (a M a ) and its inverse b = (6°^) are 
real, any matrix A^ that is hermitizing for the 7^ a, s is also hermitizing for 
the 7 M 's, and conversely. Using this fact, it is straightforward to modify the 
proof of Theorem 6 in Ref. 0] so that it applies to DFW. 

Thus, at each point X in the spacetime V, for both TRD and DFW 
theories, we have hermitizing matrices A(X) and B(X) for the sets (^(X)) 
and (ct M (X)), respectively. For TRD, the components A^ of matrix A build 
a covariant second-order tensor [I]. For DFW, we may choose A = A^. 
However, since A is uniquely defined only up to a real factor, we might also 
multiply by a positive real scalar field A: 

A = X{X)AK (16) 

Because the "flat" matrices 7^ of Eq. (fT5l) . hence also A\ are constant 
matrices in the DFW theory, it follows from ffT6|) for each point X 6 V 
that A(X) is a matrix with scalar components in the DFW theory. At each 
spacetime point, the matrix 

B = A-f° (17) 

is positive definite for both the TRD and DFW theories by the Theorem 
above. Note that, by ([I]) and ([TBI) , we get conversely 

A = Ba°. (18) 



Using Eq. ffl2l) . we define a Hermitian product between complex four- 
vectors (or quadruplets of complex scalars) u, v by setting 

{u, v) = A pv u p *v u = u^Av, (19) 

and in the same way, we define a positive-definite Hermitian product, by 
setting 

(u:v) = B pv u p *v v = u ] Bv. (20) 

In the TRD theory, the four- vectors u, v 6 TcV v, where TcVx is the com- 
plexification of the tangent space TVx at X £ V. u It results from ([12]) [resp. 



In the TRD theory, Eqs. (Tl9l) 9 and ([20|) 9 exploit an abuse of notation that views A 



from (1141) ] that each of the 7 M [resp. a M ] matrices is a Hermitian operator for 
the product ([12]) [resp. (120]) ]. that is, 

(Yu,v) = (u,rv), /i = 0,...,3, (21) 



(pfu'.v) = (u: afv), /i = 0,...,3. (22) 

Finally, we note that, due to Eqs. ([17]) and (ITS]) , we have 

(u, f ) = (a°w :«) = («: a%) (23) 



and 



(u : v) = (7° u , «) = (u, 7%). (24) 



3 Current conservation 

3.1 Definition of the current 

Obviously, the definition of the current should involve the Dirac matrices. 
There are an infinity of sets of fields (7^) that satisfy the anticommutation 
relation (jTJ) in the given curved spacetime (V, g^ v ). In the case of a flat space- 
time, it is natural to assume that the matrices 7^ are constant in Cartesian 
coordinates (such that the metric has the standard form, g^ = r/^). U In 
the latter case, the current is unambiguously defined as 

j» = ( 7 ^, vo = a pv (yx r*r, (25) 

or equivalently [using (TT2]) ]: 

j^^yt^^^t^ (B" = ^y). (26) 



and -B in two different ways: first as tensors and then as matrices representing Hermitian 
forms [4]. Compare these equations with the similar equation using the spacetime metric 
G viewed first as a tensor and then as a matrix representing a quadratic form. In this 
case we would have g pv u p v v — u T Gv for u,v € TVx- This is harmless as long as the 
covariance is being carefully checked. 

5 We note in passing that, for TRD, this is equivalent to say that the 7 M 's are covariantly 
constant: D a ^ — 0, or explicitly D n ^ p = 0. Whereas, for DFW, the derivatives D a j^ 
are always zero |12) : contrary to the Levi-Civita connection, the spin connection does not 
recognize anything special in the case of constant Dirac matrices in the flat situation. 
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Indeed, this definition coincides with the usual one [TQl EH] for the standard 
set (7^) of "flat" Dirac matrices — for which set A = 7^° turns out to be 
a hermitizing matrix. The case with the standard set has been thoroughly 
investigated in the literature, in particular the current ( 126]) [with A = 7" ] is 
then derived from the Dirac Lagrangian, so that there is no ambiguity. In 
addition, it turns out [I] that the current (|2"6"]) is actually independent on the 
choice of the Dirac matrices: if one changes one set (7^) for another one (7^) 
[satisfying the same anticommutation relation ([[]) as does (7^)], then the 
second set can be obtained from the first one by a similarity transformation: 

3SGGL(4,C): ^ = SYS-\ /i = 0,...,3, (27) 

for which the solutions of the flat-spacetime Dirac equation exchange by 

i) = Sip. (28) 

The hermitizing matrix is transformed thus: 

A = (S-^AS- 1 = (S^)- l AS-\ (29) 

and this leads indeed to the invariant relation 

p = ^^Aij; = J*. (30) 

The definition (126]) of the current is thus the right one in the flat case, and it is 
generally-covariant, the current being indeed a four-vector, for TRD and for 
DFW as well — as it results immediately from the transformation behaviours 
of its ingredients 7 M , ip, and A in the two theories. Therefore, we assume 
( 12"6"|) [or ( 12 5p ] as the definition of the current in the general case of a curved 
spacetime. We note that the invariance (130]) of the current after a similarity 
( l27]) -( ]29]) remains in force in that general case, even if matrix S depends on 
the event X. We check from (126]) . (|T7|) and the positive definiteness of matrix 
B that the probability density is J° = ip^Bip > 0, with J° > if ip 7^ 0, as 
it must be. 

3.2 Characteristic condition for current conservation 

From the current definition ([26]) and Leibniz' rule, one gets trivially 

D^ = (D^B^ + ^(D^)iP + ^B»D^, (31) 
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provided the derivative D p B^ is well defined. Note that, for TRD, B^ = 
(A'j^vp = A va r y p l(T is a spacetime tensor. U D p B^ is then the matrix with 
components D^Btf , the latter being defined in the standard way from the 
relevant affine connections, thus 

DpBS p = d,B» p + {^} BZ p - K„B% P - A^B^. (32) 

For DFW, 7^ is a "spinor vector" [Eq. f ill)) ] . whose spin covariant derivative 
is defined as [121 [35] 



D.Y^d^ + UY + iT^} (33) 

(where [M, N] = MN — NM), and this is known to be zero: 

Dtf = (DFW); (34) 

while the derivative of a "spinor scalar" like the hermitizing matrix A is 
defined to be 

D,A = d,A - AT, - rt A, (35) 

ensuring that D^(A^) = (D^Ay . Since the spin matrices have the form [12] 

T M = c Xuil s Xu (DFW) (36) 

with real coefficients c\ ull , and where s Xl/ = | (j x, y l/ — 7^7 A ), it follows from 
( 135]) and the hermitizing character ffl2]) of A that 

D M A-^A = (DFW). (37) 



6 More exactly, the first contravariant index [i in 7^ CT corresponds to tensor components 
of the (real) tangent space, as indicated by the (real) anticommutation relation ([T]). In 
contrast, the second contravariant index, c, as well as the covariant index, p, correspond 
to tensor components of the complex tangent space, as apparent from Eq. (|25p. In 
other words, the tensor 7^ dr at each spacetime point X is an element of the tensor space 
TVx <S> T c Vx (g) T£Vx, where T c Vx and T£Vx are the complexifications of the real 
tangent space TVx and its dual T°Vx- Similarly, from (fT9|) . the tensor A va belongs to 
TqVx <8> T^-Vx, so that the contracted product B^ p — A va ^ a makes sense and belongs 
to TVx ® T^Vx <g> T^Vx- Now, for different types of indices, we may use different 
connections. For "complex" indices, we are using the connection AjL, Eq. (j3|), which 
differs from the metric connection {p„} for TRD-2; but, for "real" indices, we shall use 
the metric connection for both TRD-1 and TRD-2. Leibniz' rule still applies. 
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Moreover, for DFW, we may choose A = 1 in Eq. (I16p . i.e., A = A\ so that 
we get 

D^A = (DFW, A = A*), (38) 

and from ( J34p . by using Leibniz' rule: 

D^ = Q (DFW, A = A*). (39) 

Note that, if one uses also matrices T^ for TRD, as defined from the con- 
nection coefficients [Eq. (JTJ)], then the definitions ( 133]) and ( 1331) . and in fact 
all definitions of covariant derivatives used in DFW theory [35J, also apply 
to TRD — but, of course, not in general the results (134"|) and (1391) . which are 
specific to DFW theory. 

Since A and B^ are hermitian matrices, the first term on the r.h.s. of 
(ED is 

(D^B^ = {B»D^)H = [A(YD^)]^ = {^D^Aijj. (40) 

Therefore, if the (relevant) Dirac equation (J2J) is satisfied, then the two ex- 
treme terms on the r.h.s. of ( l31~j) cancel one another: 

(^D^Aip + ^Ai'fD^i/)) = imiftAij) + tp ] A(-irmp) = 0, (41) 

whence from (13"T|) : 

D^ = ij\D,B^. (42) 

Our definition of the covariant derivative uses the metric connection for both 
TRD-1 and TRD-2 when "real" indices are concerned (see Footnote 6) — as 
is the case for the index /i in the foregoing equation. That is, for TRD, as 
well as for DFW, the connection acting on the probability current J M is the 
Levi-Civita connection. Hence, it follows that, for TRD as well as for DFW, 
the condition for conservation of the probability current is the same as in a 
Riemannian spacetime, namely D^J^ = 0. Hence we can state the following 
result: 



Theorem 1. Consider the general Dirac equation $£), thus either DFW or 
any of the two TRD equations. In order that any ip solution of |J|) satisfy 
the current conservation 

DpJ" = 0, (43) 
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it is necessary and sufficient that 

D^B* = (£" = Af). (44) 

The reasoning which was used to get Theorem 1, hence also this theorem 
itself, extend immediately to the transition probability currents 

K^{^,^')=^B^'. (45) 

We note also that Eq. (1391) gives to Theorem 1 the following 



Corollary 1. For DFW theory, the hermitizing matrix field can be imposed 
to be A = A\ with A^ a constant hermitizing matrix for the "flat" constant 
Dirac matrices rf a of Eq. (T73J). Then the current conservation applies to 
any solution of the DFW equation. 

The current conservation usually stated for DFW theory {e.g. [i~2 | l24" t l29]) 
applies only to the case where the "flat" Dirac matrices 7^ a of Eq. ffT5|) are 
the standard choices of Dirac matrices, which we shall denote generically 
7"°, and for which 7" turns out to be hermitizing. {These choices of Dirac 
matrices are related by similarity transformations that are unitary [22], and 
not by general similarity transformations f[27|) that leave the anticommutation 
formula ([1]) invariant.} This is a particular case of Corollary 1. 

3.3 Modified Dirac equation with conserved current 

Due to the r.h.s. of (I42p . the current is in general not conserved for solu- 
tions of the Dirac equation (T5]) (the general one, i.e., DFW or TRD as well, 
although for DFW the natural choice A = A^ does ensure the current con- 
servation). However, we can modify this equation to the following one: 

7"Z>^ = -irmp - ^A-\D^)iP, (46) 

so as to conserve the transition probability current (T4"5"|) . Indeed we have the 
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Theorem 2. For any pair (ip,ip') of solutions of the modified Dirac equa- 
tion (Jfi), the transition probability current (J5) is conserved. 



Proof. Similarly to Eqs. (13TJ) and (I40p , we get from the definition fj4"5]) : 
D^ = (-fD^Atf + ^(i^B'V + ^A'fD^'. (47) 

Then, if the modified Dirac equation ( I46p is satisfied by ip and by if)' as well, 
the contributions coming from —imip and —irmp' to the two extreme terms 
on the r.h.s. of ( l4"T|) cancel one another, as in Eq. ( |4"T|) . Thus we are left 
with 

D^K" = (C0Aip' + ^{D^B")^' + ip^ACip', (48) 

where 

C = - l -A-\D^). (49) 

But since A = A\ we have {A^ 1 )^ = A" 1 , hence [noting that (D^B^y = 

{C^ = tftf = - l -^{D,B»)A-\ (50) 

so that the r.h.s. of (l48l) vanishes, as claimed by lTheorem 21 Q.E.D. 



Note that the modified Dirac equation ( 1461) coincides with the normal 
one ([2]) for all ip, iff the condition for current conservation ( l4"4"|) is satisfied. 
Therefore, Eq. (j4^1) is really the adequate modification of (J2J) to get the cur- 
rent conserved in the general case. However, it is the normal Dirac equation, 
not the modified one, that (in the TRD case) has been derived from the 
classical-quantum correspondence [3]. Moreover, the TRD-1 version [based 
on the Levi-Civita connection (J3J)] of the normal equation obeys the equiva- 
lence principle in the precise sense of the IIntroduction( whereas it is not the 
case for the corresponding version of Eq. (|4U|) , because the vanishing of the 
connection (jlj) does not imply the validity of the condition (I44j) . Thus, one 
may feel that the physically relevant equation remains the normal one (J2J). 
Since the current conservation is very important, this option means that not 
all possible fields 7^, A are physically admissible, but merely the ones which, 
in addition to the anticommutation relation (QJ), satisfy condition fi44\ ). Such 
systems will be called admissible. This, after all, is just an extension to the 
general case of the statement made for the flat case, that any relevant field 
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7 M (and hence also the field A) has to be constant in Cartesian coordinates: 
if one selects the gamma field at random, the condition (j44p and the current 
conservation do not generally apply to the solutions of ([2]) even in a flat 
spacetime — except for DFW. 

3.4 Similarity transformations under which the Dirac 
equation is covariant 

Consider a local similarity transformation of the coefficient fields: starting 
with a "fiduciary" set of fields, 7^, A, B^ = A^, let us apply an event- varying 
similarity transformation S(X) [we exchange S for S^ 1 w.r.t. f )27|) -( l29|) . for 
convenience] : 

f = £- VS, (51) 

A = S^AS, B» = AY = S^B^S. (52) 

For DFW, the similarity transformations must be restricted to the form S(L) 
with L i—)- ±S(L) the spinor representation, see Eq. ( 1105(1 . As is well known, 
the DFW equation is covariant under all such similarities [12] . For TRD, 
we may ask whether the general Dirac equation — either the normal one (J2J) 
or the modified one (|4~6I) — is covariant under a such transformation, if one 
simultaneously imposes that the wave function must transform naturally as 

# = S^ip. (53) 

We examine three questions: i) ip obeying the modified equation, when does 
ip obey the normal one? ii) ip obeying the normal equation, when does if) 
obey it also? and iii) ip obeying the modified equation, when does ?/> obey 
it also? The first question occurs most naturally after the discussion at the 
end of subsect. 13.31 The answer to it is given by 



Theorem 3. In order that ijj obey (TJj) each time that ip obeys (pfy , it is 
necessary and sufficient that 

Z = B^D^S + \{D^)S = 0. (54) 
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Moreover, if Z = 0, the transformed fields satisfy condition fi44\ ), hence the 
current conservation is valid, after the transformation, with the normal Dirac 
equation (TJj]. 

Proof. Entering (153]) into (146]) . we get successively 

= S[(S- 1 YS)D fl + tm]4, + Y(D ll S)4> + ^A- 1 (D^)S^ (55) 
= SVip + A^Z^, (56) 

the latter because A~ X B^ = 7 M , and where T> is the normal Dirac operator 
with the transformed fields j^, the normal Dirac operator being T> = 7^.0^ + 
im. This proves the first part of Theorem 3. It remains to check that 

Z = => D ft B' t = 0. (57) 

Entering into ( 1571) 9 the definition ( 152]) of B^ as a transformed quantity, yields 
S^B^D^) + (D^)B^S = -S^D^B^S. (58) 

Thus, ( 157|) 9 is equivalent to 

Y s = - l -S\D^)S, (59) 

where 

Y = S^B^D^S, (60) 

and where Y s = ^(Y + Y"l") denotes the Hermitian part of Y. But the r.h.s. 
of (1591) is a Hermitian matrix. Hence the equation Z = 0, which is just 

Y = - l -S\D^)S, (61) 

is equivalent to the conjunction of (159"]) and the vanishing of the antihermitian 
part of Y, namely 

S^B^D^S - (D^S^B^S = 0. (62) 

In particular, Z = implies (159]) . and thus implies (ISTI) ?. This completes the 
proof. 
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Thus, coming back from the modified Dirac equation to the normal one 
leads necessarily to a normal equation for which the current is conserved. 
Note that Eq. ([5"lj) is a first-order linear system of 16 independent PDE's 
for the 16 independent unknowns of the similarity matrix S(X). So that, 
starting from fields 7 M , A "selected at random," solving this system (which 
should be possible with suitable boundary conditions) allows us to go to ad- 
missible fields 7 M , A. The meaning of the additional condition ( 15"2"|) will be 
soon clarified by ITheorem 41 For now, recall that, if we start from a system 
of fields 7 M , A that already satisfies condition (j4"4"|) . then the "modified" Dirac 
equation actually coincides with the normal one. Hence, ITheorem 31 allows 



us to immediately answer question ii) put at the beginning of this subsection: 



Corollary 2. Let us start from a set of fields 7^,^4,5^ = A^ that does 
satisfy the condition ft44\ ) for current conservation. In order that ip still obey 
(TJ|) each time that ip obeys it already, it is necessary and sufficient that 

Z = B^D^S = 0. (63) 

Moreover, condition [J4\ ) is then preserved by the similarity transformation. 

Note that Eq. (1631) is the particular case D^B^ = in Eq. (154")) . and 
is thus also a first-order linear system of 16 independent PDE's for the 16 
independent unknowns of the similarity matrix S(X). 

Finally, the answer to question iii) is given by 



Theorem 4. In order that ip still obey the modified equation fijffy , each time 
that ip obeys it already, it is necessary and sufficient that Eq. [6^1 be satisfied. 

Proof. Let us write that ip obeys Eq. (J4U1) . and multiply this by S: 



= S 



tfDrf + %m^ + l -A-\D^)i, 



(64) 



On the other hand, the fact ip obeys Eq. (J4U1) is equivalent to Eq. ( )55|) . thus 
to 

= S [fDr + im]iP + ^(D^ + ^A-\D^)Si>. (65) 



Hence, ip obeys (jUJ) at each time that ip obeys it already, iff: 

V^ ^(D,S)ij + 1 -A-\D^ 1 )S^j = ^SA-\D,B^. (66) 

Inserting the expression ( |52i) of A and B^, and thus computing D^B^, yields 
SA-\D^) = A-\D^)S + Y(D^S) + A-\tf)- l (D^)B»S. (67) 
Thus, the characteristic condition ( 166]) rewrites as 

Y(D^S) = A-\S^-\D^)B»S, (68) 

or (remembering that A^ = B^): 

&B»(DrS) = (D^)B^S, (69) 

which is precisely Eq. ( 1621) . Q.E.D. 



Let us summarize. If a similarity transformation takes the modified Dirac 
equation ( 146]) to the normal one ([2]), then it leads to a normal equation with 
admissible fields 7^, A, i.e., ones for which, in addition to the anticommmu- 
tation relation (CQ), the condition ( 1441) for current conservation is satisfied 
( [Theorem 31) . And if, starting with admissible fields 7^, A, the normal Dirac 



equation is covariant under the transformation, then necessarily the trans- 



formed fields are admissible also (Corollary 2) This provides the justifica- 



tion for the restriction to admissible fields. Finally, the condition ( 1541) of 
ITheorem 31 that allows one to come back from the modified to the normal 
Dirac equation, turns out to be equivalent to the conjunction of the condition 
( l4~4"j) for current conservation, and of condition (1621 . The latter is just the 
one ensuring the covariance of the modified Dirac equation (ITheorem 4j) . 

4 Hermiticity of the Hamiltonian 

4.1 The Hamiltonian operator and its frame depen- 
dence 

The general Dirac Hamiltonian H is obtained by multiplying the general 
Dirac equation (|2]) by 7 on the left, using the anticommutation formula ([I]). 
This puts the Dirac equation into Schrodinger form: 

^ = H ^> (t = x% (70) 
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with 

H = ma° -ioPDj-i{D -d ), (71) 

and where the a M 's are given by Eq. ({TBI . One sees from the latter equation 
that, for TRD, a^f = (a") p u is not a general tensor, in contrast with 7^ p = 
{l^Yv — an< i that, however, it does behave as a spacetime tensor for purely 
spatial transformation of coordinates: 

X* = x°, z"' = f((x k )). (72) 

This is natural. Indeed, the rewriting of any linear wave equation in the 
Schrodinger form (1701) . in which the (linear) Hamiltonian operator H has to 
contain no time derivative, is based on a splitting of spacetime into space 
and time. Thus, a priori, the operators H and H' corresponding to different 
spacetime coordinate systems should be different, in general. This can be 
checked from the transformation of the Schrodinger equation f[7U|) . Consider 
first the case that the wave function is a scalar, or transforms as a scalar 
(which is relevant to DFW). Then Eq. ( 170]) transforms simply as 

.dip' .dipdx" dt ,dipdx j 

H-w = z— — = i— — = — — Hip + i— — ——. 73 

Y dt' dx» dt' dt' Y dxi dt' v ; 

In order that the Hamiltonians H and H' be equivalent operators, it must 
be that: i) dx 3 /dt' = 0, ii) Hip transforms as a scalar: (}i'ip')((x' 11 )) = 
(H.ip)((x v )), and iii) t! = t, that is, x'° = x°. l\ The three foregoing condi- 
tions impose that we restrict the allowed transformations to be purely spatial 
coordinate changes (172|) . Now, for DFW, H^ as obtained from Eq. (171 P is 
indeed a spatial scalar [2]. Next, consider the case that the wave function is 
a four- vector (which is relevant to TRD). We have then 

,„, mu .d<ip'» . d (dx'» IV \ dx'" f.dip v dxP\ 

the latter for linear coordinate changes. Thus, in order that the new Hamil- 
tonian H' that appears in Eq. (1711) be an equivalent operator to H, the 
expression above should now transform as a four- vector. This needs that 



7 Since we set t = x° (and t' = x'°), the "true" time is rather T = x°/c. Note that 
x° = cT is invariant under a change J" = aT. In other words, time scale changes are 
(fortunately) allowed. 
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or 



That is, again one must restrict oneself to spatial changes (1T21) . Since this 
restriction applies already to linear coordinate changes, it must be imposed 
to all coordinate changes. One verifies that, under all spatial coordinate 
changes (I72"|) . H'ijj', thus defined as the transformation of the l.h.s. of the 
Schrodinger equation (|70|) . is indeed the four- vector transformation of H0. 
One then easily checks that Hip, as defined instead from the explicit expres- 
sion ( ITT]) of the Hamiltonian, does transform as a four-vector under spatial 
changes ( 172|) for TRD. Thus, for DFW and TRD as well, in order that the 
Hamiltonians H and H' before and after a coordinate change be equivalent 
operators, the coordinate change must be a spatial change (I72"|) — then, both 
sides of the Schrodinger equation (170)) behave as a scalar for DFW, and as 
a four-vector for TRD. This is consistent with the fact that, both ip and 
Hip being wave functions in the same Hilbert space, both ip and Hip must 
transform the same way under "allowed" coordinate transformations. 

Thus, the Hamiltonian operator associated with a given wave equation 
depends on the reference frame F which is considered — the latter being un- 
derstood here as an equivalence class of local coordinate systems (charts) on 
the spacetime V, modulo the purely spatial transformations (172]) . g Since 
the time coordinate x° is fixed, this is a more restrictive definition of a ref- 
erence frame than in relativistic gravity, where changes x'° = /((x M )) are 
allowed [13] [23] . However, one may trace back this restriction to that ef- 
fected by mechanics itself: e.g. in an inertial frame in a flat spacetime, the 
inertial time (synchronized according to the Poincare-Einstein procedure) is 
naturally distinguished; accordingly, the quantum Hamiltonian will generally 
change if one selects another time coordinate, see Eq. (173"]) . As shown in Ref. 
[6] , the data of a reference frame F determines a three-dimensional "space" 
manifold M, which is the set of the world lines of the observers bound to 



8 This notion of a reference frame is formalized in Ref. [6 , together with the notion 
of the associated space manifold M, which is time-independent. This formalization needs 
that one restricts oneself to an open domain U in V, such that there is at least one chart 
X ■ U — > R 4 , thus that covers U. Since a chart is, from its definition, a difieomorphism of 
its domain onto its range in the arithmetic space R 4 ~ R x R 3 , this means that we assume 
a local 1x3 decomposition of spacetime — as is indeed commonly assumed in works on 
quantum theory in a curved space-time [TS]. This decomposition writes X i-> (x°,x) 
with x £ M: see Sect. 4 in Ref. [§], point i). 
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F — i.e., whose spatial coordinates x J do not depend on the time x°, in any 
chart of the class F. The charts of the class F, also called charts adapted 
to F, provide an atlas of M: the coordinates of the running element of M in 
such a chart are just the constant spatial coordinates w* of the corresponding 
world line. 

Equation (I7T|) reveals an important property of the TRD-2 connection 
d5J) in the preferred frame E. Namely, using the definition (jSJ) to set Dq = do 
in Eq. (1711) shows that the Hamiltonian operator for TRD-2, at each time t, 
depends only on the choice of gamma matrices and the spatial geometry at 
time t. 



4.2 The scalar product 

Ideally, we would like to define in a natural way a Hilbert space scalar prod- 
uct for wave functions, and to find that "the Hamiltonian H is Hermitian 
for this scalar product." We now know that this property, like H itself, will 
likely depend on the reference frame, so we select one, F, with the associated 
"space" manifold M. Note that H, because it does not involve time deriva- 
tives, primarily operates on spatial wave functions if) = ip(x) with x G M. 
However, in general, H and the metric g^ v do depend on the time t = x°. As 
noted in Ref. [2J, the scalar product for the Dirac equation ([2]) has necessarily 
the following general form: 

(V \<p) = [ ty{x).<p{x)) dV(x), (77) 

where (u.v) is a Hermitian product defined for "arrays u and v of four complex 
numbers" (which in fact are either complex four-vectors or quadruplets of 
complex scalars, depending on whether TRD or DFW theory is considered), 
and where dV(a;) at each time t is an arbitrary volume element defined on 
M. The latter has the form (at time t) 

dV(ar) = a(t, x) y/-g(t,x) d 3 x, g = det(^) (78) 

where a(t,x) is any spatial scalar field [i.e., scalar under transformations 
(172]) ]. Note then using formula f!72|) that the integral defining the scalar 
product (I77|) is invariant under all coordinate transformations of M. Also 
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note that the volume element a(t, x) yj—g(t,x) d 3 x is the most general pos- 
sible on M. 

For a flat spacetime, in any Cartesian coordinate system (x^), the follow- 
ing scalar product has been identified: 

(V> \\<p) = / (V'W : p(x)) d 3 x, (79) 

J space 

where [u : v) is the positive-definite product ( ]2"0)) with the corresponding 
constant hermitizing matrix B. It is fully satisfying in that case, since the 
Hamiltonian is always Hermitian for that product [3]. Hence, the sought 
after scalar product ( 177|) must coincide with (179)) in a flat region of U, when 
it is endowed with Cartesian coordinates. More generally, for any given 
event X e U, and for wave functions of a priori bounded variation, vanishing 
outside a small neighborhood of X, the metric and the matrix B may be 
considered constant, so that one should be able to approximate the exact 
product (I77|) by the product (1791) . rewritten in a covariant form. Therefore, 
we must take (i/}(x).<p(x)) = (ip(x) : tp(x)) in (177)) . and we are left with 
the mere choice of the volume measure V on M, thus with the choice of the 
scalar o. An obvious possible choice, indeed the standard choice, is a = 1. 
So that the general form of the possible scalar product may be written, both 
in intrinsic form and in a chart x on the space M, as 

(ip\(p)= (ip : <p) dV = / ^Aj°cp o v^ d 3 x. (80) 

JM ^x(M) 

For DFW theory, recall that the matrices 7 M , in particular 7 , are defined 
by Eq. (JT5J) , in which the matrix (a M a ) must be invertible and satisfy 

a\a»^=g^. (81) 

If we make the standard choice a = 1, and take for "flat" matrices the 
standard Dirac matrices 7'^, so that we may choose A = 7^° as previously 
noted, then, in that particular case, the scalar product (180)) is the one stated 
by Leclerc [23]. If, in addition, we have a • = 0, we get a° = \J g m from 
USD, hence 

v ^A 7 ° = v /Z ^l4. (82) 
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Now assume that g^j = (j = 1,2,3), as is in particular the case for a 
static metric in adapted coordinates. Then we can certainly get a°j = 0, and 
moreover we have g 00 = 1/goo, whence 

^]Aj = VhU, (83) 

where h = det(hjk) is the determinant of the metric hj^ = —gjk {j — 1 5 2, 3) 
on the space M, induced by the spacetime metric g^ u . Thus, in the case that 
9oj — (j = 1, 2, 3), and using furthermore the standard Dirac matrices, the 
definition (JED} coincides with that found "natural" in Ref. [2] [Eqs. (21)-(23) 
there] for a static metric in DFW theory. 

4.3 Conditions for hermiticity and isometric evolution 

In quantum mechanics, the Hamiltonian operator H(t) is an operator-valued 
function of time t, that represents the energy observable at time t. There are 
three axioms in quantum mechanics that can be straightforwardly adapted 
to a curved spacetime as follows: 



Axiom (A). The Hilbert space scalar product (8(fy of any time-independent 
wave functions ip and ip defined on the space manifold M is time independent. 



Axiom (B). For each time t, the Hamiltonian H is a Hermitian operator 
with respect to the scalar product ( fffOI) . 



Axiom (C). The solutions of the Dirac equation $7D ) have an isometric 
evolution with respect to the scalar product (E2P- 

It will be shown in this subsection that Axioms (A) and (B) together 
uniquely determine the Hilbert space scalar product, up to an inconsequen- 
tial constant. 

Note that these axioms are somewhat weaker than the self-adjointness 
and unitary evolution required in quantum mechanics. Axiom (B) is a pre- 
requisite for the energy operator H(t) to be self-adjoint for each time t. Axiom 
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(C), which guarantees that the Hilbert space norm is preserved in time, is a 
prerequisite for ensuring that the solutions of the Dirac equation (ITU]) have 
a unitary evolution with respect to the scalar product ( 180]) . 

It can be shown that any pair of these axioms implies the third. In 
particular, we note that if Axiom (A) is valid, then Axioms (B) and (C) are 
equivalent, as they are in ordinary quantum mechanics. This observation 
follows immediately from differentiating Eq. (ISO]) with respect to time to 
obtain 

d (ip \<p)= /V d<>[oJ=gAf\ ip d 3 x + ((d ip) \<f) + (ip\ (d <p)) (84) 



M 



and then substituting for if and if either time-independent wave functions or 
time-dependent solutions of the Dirac equation ( 170]) . In the former case, from 
Axiom (A) we have do(ip \ f) = 0, and since if and if are time- independent 
we get 

= I ^ dolay/^A^ } if d 3 x + 0. (85) 

Jm 

Since if and f can be chosen as arbitrary smooth, 4-component complex 

wave functions with compact supports in M, the vanishing of that integral 

implies that 

d [ay^A 7 °]=0. (86) 



Then, in the latter case, for two solutions ip and if of the Dirac equation (I70|) , 
we may substitute Eq. (186]) and H = ido in Eq. (184]) . obtaining 

d Q (ip\if)=i[{Kip\if)-(ip\Hif)\. (87) 



Thus, assuming Axiom (A) , the hermiticity of the Hamiltonian is equivalent 
to the isometricity of the evolution, generalizing a well known result in quan- 
tum mechanics. 

The hermiticity condition for the Hamiltonian is, by definition, 

(HV \if) = (V I Hp) (88) 

for all if and if in the domain of H(t), denoted Dom(H), which we assume 
to be independent of time. Recall that the wave functions if and if in this 
definition are time-independent. Because the sesquilinear form (Hip \ if) 
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is determined by the corresponding quadratic form Q(ip) = (Hip | ip), an 
equivalent condition to (I8"8"j) is actually 

Vip G Dom(H) (Hip \ip) = (ip\ Hip) (for each time t). (89) 



In order to find when this does occur with the Hamiltonian (ITT]) , we use the 
expression of the global scalar product (|80|) as an integral of the Hermitian 
product (u : v ), and we use the relation between the two Hermitian products 
(u : v) and (u,v) in Eq. (12~4"|) . exploiting their Hermitian property to obtain: 

(ip : HV>) - (Hip :ip) = (ip : -ia j Djip) + (ia j Djip : ^) 

+ (</> : -i(A) - d )ip) + (i(D - d )iP : iP) 

= (ip, -i^DjiP) + (ijWjiP, iP) + (ip, -ij°D ip) 
+ (tf°D iP,iP) + (1P : id iP) - (td iP : iP) 

= (iP,-irD^) + (tYD^iP,iP) 

+ (iP : id iP) - (id ip : ip). (90) 

Therefore, if the wave function ip is a time-dependent one that obeys Dirac 
equation in the Schrodinger form f lTU]) . we get 

= (ip, -i-fD^) + (ij^D^iP, iP), (91) 

which is indeed an immediate consequence of the Dirac equation in the initial 
form ([2]). Thus we are left with = in that case. 

However, if the wave function ip is a time-independent one instead, as is 
normal at the stage of checking the hermiticity of the Hamiltonian, what we 
get is more interesting: 

(iP:HiP)-(Hip-.iP) = (iP,-i^D^) + (iD^,^iP) 

-i[^A^D fi iP + (D^A^iP}. (92) 

Now, recall that the field 7 M , A may not be selected freely, but instead should 
satisfy the admissibility condition (jUJ) [and relation (0Q)]. If this is the case, 
we find thus [cf. Eq. (12%])]: 

(iP : Hip) - (Hip : ip) = -iD^A^ip) = -iD^. (93) 

As noted after Eq. ( l4"2"j) . for TRD, as well as for DFW, the connection acting 
on the probability current J M is the Levi-Civita connection, hence 

v^ IV" = d, {y/=jj J») . (94) 
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We use the form flSDJ) of the scalar product. With this scalar product, by 
integrating (193 j) over the space manifold, we get (setting the boundary term 
equal to zero by assuming that the functions ip G D decrease quickly enough 
at spatial infinity): u 



z[(V> | Kip) - {Kip | ip)} 



a ^f^jD^d 3 * 



(95) 



R :i 



[ ad,(^gJ») d 3 x 
/ [a d (y^ J°) - (d j( r) 

JR3 



-p J j ] d 3 x + / dj {(jV^9J J ) d 3 

JR3 



R 5 



H^) - {Kip 



R 5 



[o" <9 (v^^WV) - {djff) s/^gip^A^ip] d 3 x 
V> f [a d {yf^Af) - (V) 4^9 V] ^d 3 x. 



Using again the fact that a sesquilinear form is determined by the correspond- 
ing quadratic form, we find thus that the hermiticity of the Hamiltonian is 
equivalent to ask that, for all ip and <p in Dom(H), 



ip ] Nip d 3 x 



0. 



N 



a 



R3 



do (v^7°) - (<V) V=g An*. (97) 



In the same way as after Eq. ( |85|) . the matrix A" must thus vanish for every 
x G R 3 . Hence, the characteristic condition of hermiticity for the general 
form (ED]) of the scalar product is 



o" <% {^gA^) - (djcr) y/=g Aj j = for every x G R 3 



(98) 



This result opens the possibility that, in a given metric and with a given ad- 
missible set of fields 7^, A, one might get the Hamiltonian Hermitian by an 
appropriate choice of the scalar field a, that determines the scalar product. 
However, combining Eq. f )98|) with Eq. f )86|) shows that {y/— g A^) d^a = 0. 
That is, since the matrix y/—g A is invertible, ^d^a = . Then, since the 
gamma matrices are independent, d^a = 0, proving that the scalar field o is 



indeed constant, assuming Axioms (A) and (B) Without loss of generality 



we may henceforth set a = 1. Thus we can state the 



9 From now on, for definiteness, we shall assume that x(M) = R 3 (for the chart x€F 
that we consider), and hence that the "space" manifold M is diffeomorphic to R 3 . 



(96) 
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Theorem 5. Axioms (A) and (B) uniquely fix the scalar product to be 



(iff\<p)= [ tfAfip V=9 d 3 x. (99) 

Moreover, rewriting (1981) with a = 1 yields the following result: 



Theorem 6. Assume that the coefficient fields 7 M , A satisfy the two admis- 
sibility conditions (Q]j and fi44\ ). In order that the Dirac Hamiltonian (71) be 



Hermitian (at time t) for the scalar product ( fgPj) . it is necessary and sufficient 
that 

iV(x) = <9 M(x) = for every x G R 3 , M = y/=g~Aj°. (100) 



The matrices in Eq. fllOOl) depend a priori on the chart, say \. However, 
as we saw, the Hamiltonian (I7T1) and the scalar product flHUl) depend only on 
the reference frame F, i.e., on the equivalence class to which the chart be- 
longs, modulo a relation of the kind "x^x' iff X an d x' exchange by a purely 
spatial transformation (!72|) " [6]. The transformation of the matrices in Eq. 
( HOOP is consistent with this: it is easy to check that the validity or invalidity 
of the condition f llOOp is invariant under any spatial coordinate change (!72|) . 
thus it depends only on the reference frame F. 



4.4 Effect of an admissible change of the coefficient 
fields 

In a given spacetime (V, g^), there are infinitely many coefficient fields 
(7 M , A) that satisfy the two admissibility conditions (0Q) and ( |44|) . The ques- 
tion thus arises, whether or not (in a given reference frame F) the hermiticity 
condition f llOOp is preserved by a change of the admissible coefficient fields — 
such changes will be called admissible changes of coefficient fields. We have 
the definite answer to this question only for DFW. For both DFW and TRD, 
we know that, if a change of the fields 7 M respects the anticommutation re- 
lation (IT]), then this is a similarity transformation ( 15~TT) . and that then the 
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matrix A and the field B^ change according to (1521) . Therefore, the matrix 
M of lTheorem 61 changes like A and B^\ 

M = y/=gAf = S j MS. (101) 

As an immediate consequence of this and ITheorem 61 we have the 



Corollary 3. Let (7^, A) be a set of admissible coefficient fields. Let us fix 
the reference frame, so that the Hamiltonian operator and the scalar product 
are determined by the fields (7^, A). Let S = S(x) be a space- dependent ma- 
trix such that condition frffl) is still valid after the transformation / f37]) -/ f5^) 
of the coefficient fields. Then, the initial Hamiltonian is Hermitian iff the 
transformed Hamiltonian is Hermitian. 

In short: the hermiticity is preserved by admissible changes that do not 
depend on time. However, admissible changes may well depend on time. 
For DFW, the admissible changes are restricted to local similarity transfor- 
mations belonging to the spin group, S(X) G Spin(l,3). o We know from 



Corollary 1 that all such transformations, including time-dependent ones, are 



admissible for DFW, since they preserve the conditions ([15]) . OS]), as well as 



1 Any such similarity is obtained from a change of the tetrad field u a = a tl a 9 M , by a 
(proper) local Lorentz transformation L = L(X) E SO(l,3): 

M^iy^Ly^^oy,, (102) 

thus 

^ p = a» e L c p . (103) 

This allows us to define gamma matrices from the same fixed set of "flat" ones through 



these two tetrad fields [cf. Eq. (fT~5|) ] : 

7^=o" a 7 l|a , 7^ = 5^7^, A = A = A\ (104) 

Using (|103|) and (|104j) together with the characteristic property of the spinor representation 
L h- > ±S(i) (defined up to a sign), it is easy to show that we have 

7^ = a" Q S - y Q S = S^^S, S = ±S{L). (105) 

Note that L and S = ±S(L) can depend on the event X, thus in particular on the time 
t. Moreover, with the restriction of S to the spinor representation: S = ±S(L) with 
L e SO(l, 3), one may show that S 1 ' A^ S = A\ thus (fTD4]| q is compatible with (|q2]) . 
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([T]). Thus, the hermiticity condition is always (llOOp . for DFW (with the 
choice A = A*). Therefore, the hermiticity conditions before and after an 
arbitrary similarity transformation S, are respectively 

d M = for every x G R 3 (106) 

and 

d (S*MS) = for every x e R 3 . (107) 

Obviously, if the initial matrix M is independent of time, then, with a time 
dependent similarity transformation S, the transformed matrix M = S^MS 
will in general depend on time — thus contradicting OqM = 0. 

To take a sensible example, consider the very general case of an admissible 
coordinate system, i.e., a coordinate system such that goo > and such that 
matrix (gjk) is negative definite. In that case, the lemma in Appendix B of 
Ref. [1] shows that the matrix (a M a ) of Eq. f lHTj) (which is matrix M in Ref. 
[I]) may be chosen to satisfy a • = 0, so that, taking for "flat" matrices the 
standard Dirac matrices 7" M , with A = 7" , we have [Eq. 



m = V-gg^u (108) 

Thus, the hermiticity condition with that starting system, Eq. (I106p . rewrites 
as Leclerc's condition 



do(V~gg 00 ) = f o r ever y x e R 3 , (109) 

while in the transformed system, it is (I107p . or here: 

d (^-gg 00 S^S) = for every x G R 3 . (110) 

Now, if ( I109p is satisfied, and if S is such that S^S does depend on time, it 
follows that (II 101) cannot be satisfied. This proves that, for DFW, the her- 
miticity condition is not invariant by the admissible changes of the coefficient 
fields {'f,A). 

For TRD, the condition ()44l) is much more demanding than it is for DFW, 
because: i) The derivatives D^ u have no reason to vanish in general, and 
ii) also the hermitizing matrix A is not covariantly constant in general. For 
instance, in the case of a flat spacetime, the condition D^ v = means (for 
TRD) that the gamma matrices are constant in Cartesian coordinates. In 
contrast, for DFW, the condition D^ = is always satisfied, so that in 
Cartesian coordinates the gamma matrices need not be constant. 
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5 Summary and conclusion 

In this work, we studied simultaneously three versions of the Dirac equation 
in a curved spacetime: the standard or "Dirac-Fock-Weyl" (DFW) equa- 
tion, and two tentative versions, proposed recently |3j . In a given coordinate 
system, the three equations differ merely in the covariant derivative. The 
two alternative versions are based on the tensor representation of the Dirac 
(TRD) field, Eqs. (JSJ) and (jUJ). In order to define a conserved probability 
current for any one of these Dirac equations in a curved spacetime, we must 
consider a general set of Dirac gamma matrices (7^) and a hermitizing matrix 
A. We call these matrices (7^, A), which vary with each spacetime point, the 
coefficient fields of the Dirac equation. Different choices for the coefficient 
fields (7^, A) are related by local similarity transformations. For the two al- 
ternative equations based on TRD, these local similarity transformations are 
not restricted to Spin(l,3) transformations (associated with Lorentz trans- 
formations of a tetrad [12]), nor to the unitary transformations which are 
considered in Ref. [52], but instead comprise the entire group of GL(4, C) 
transformations. 

Independently of which Dirac equation is selected, the current conserva- 
tion asks for an admissibility condition to be satisfied by the coefficient fields 
(7 M ,A), as shown in ITheorem 1[ This condition restricts only the choice of 



A for DFW and is thus essentially always verified for DFW (Corollary 1) 
But it also strongly restricts the choice of the field 7^ for TRD. However, 
one may modify the gravitational Dirac equations (any of them) so that the 
current conservation is always satisfied, as proved in ITheorem 21 Moreover, 
starting from any coefficient fields (7^, A), one may transform to (7^, A) sat- 
isfying the admissibility condition, by a local similarity transformation — this 
is shown by ITheorem 31 For TRD, the local similarity transformations un- 
der which the "unmodified" or the "modified" Dirac equation is covariant 
are characterized by |Corollary 2| and ITheorem 4| respectively. 



We prove in ITheorem 51 that the Hilbert space scalar product is fixed by 
the axioms of quantum mechanics. Note in passing, that some authors choose 
a Hilbert space scalar product that differs from Eq. (1991) . and therefore their 
Hilbert space does not satisfy the axioms of quantum mechanics listed as Ax- 
ioms (A) (B) , and (C) in Section 14.31 [7]. The hermiticity of the Hamiltonian 



for the scalar product compatible with the axioms is studied in the general 
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case in lTheorem 6\ i.e., for a general coordinate system in a general curved 
spacetime. In particular, in any time-independent metric, the Hamiltonian 
is Hermitian for any admissible choice of time-independent coefficient fields 
( / y fl ,A). However, for the standard equation (DFW), the Hamiltonian in a 
given coordinate system may be Hermitian for some admissible choice of the 
coefficient fields (7 M ,v4), and non- Hermitian for another admissible choice — 
see Subsection 14.41 It means that the very existence of an isometric evolution 



of the states (Axiom (C) ) depends on an arbitrary choice of coefficient fields 
(7^, A)— for DFW. 

One may ask if the use of the TRD version of the Dirac equation leads to 
a clearer and more transparent presentation of the standard DFW formalism. 
The answer is yes. Note that the quantum mechanical treatment of the Dirac 
equation today is mostly limited to stationary spacetimes [18]. The reason 
for this, as discussed in Section HI is that in a non-stationary spacetime, the 
time evolution of free particle wave functions is generally not unitary, since 
the Dirac Hamiltonian is not generally Hermitian. By investigating a wider 
group of transformations of the Dirac equation, the possibility exists that 
one can extend quantum mechanics to non- stationary spacetimes, so that 
the time evolution of free particle wave functions is unitary. To explore this 
the wider group of local similarity transformations is a key ingredient. 
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